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ABSTRACT 


A  branch  and  bound  algorithm  for  finding  the  minimal 
cost  symmetric  assignment  is  discussed.  The  matching 
problem  in  graph  theory  and  the  Chinese  postman  puzzle 
are  all  special  cases  of  the  symmetric  assignment  problem, 
and  hence  this  algorithm  can  be  applied  to  solve  them. 


Introduct ion:  The  well  known  assignment  problem  is  to 


n  n 

minimize  7  =  T  T  c  .  .  x  .  . 

i-1  j-1  U  'J 


subject  to  Z  x..  =  1,  for  j  =  1 . n 


i  =  l 


•J 


Z  x..  =  l,for  i  =  l . n 

j»l  U 

x. .  >  o 
•J  - 
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where  n  is  a  given  positive  integer  and  C  =  (c.j)  is  a  given  n  x  n  cost 
matrix.  denotes  the  transpose  of  C  .  Any  feasible  solution  to  the  assignment 


problem  can  be  represented  by  the  matrix  X  =  (x .  j )  where  the  x„  satisfy 
conditions  (l). 


Consider  a  matrix  X  =  (x . j )  of  order  n  x  n  ,  in  which  only  one  of  the 
elements  in  each  row  and  column  is  equal  to  l,  while  all  the  remaining  elements  are 
zero.  Such  a  matrix  is  a  permutation  matrix  and  is  called  an  ass i gnment .  Every 
assignment  is  a  basic  feasible  solution  of  the  solution  set  of  (l)  and  vice  versa. 
We  use  the  letters  a,  b  to  denote  assignments. 


Occasionally  it  is  convenient  to  denote  an  assignment  by  the  set  of  its  unit 
cells,  i.e.,  the  cells  in  the  matrix  X  representing  the  assignment  which  have 
unit  entries  in  them.  All  the  other  cells  have  zero  entries  in  them  of  course.  Thus, 

a  ~  f  ( ' »  j  j ) »  ■ • • >  ( n  >  j n ) }  (2) 

where  j^,  ....  j  is  a  permutation  of  the  numbers  1,  2,  ....  n  i s  an  assignment. 
Correspondingly  we  shall  write  (i,  j)  e  a  or  (i,  j)  £  a  to  indicate  that  in  the 


/ 
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matrix  X  representing  the  assignment  a  »  the  entry  in  the  cell  (i,  j)  is  one 
or  zero  respectively. 

Let  denote  the  cost  corresponding  to  assignment  a  ,  w.r.t.  the  cost 

matrix  C  .  Hence  for  the  assignment  in  (2) 

n 

Zr  (a)  =  I  c  . 

L  r=l  rjr 

when  there  is  no  ambiguity  about  the  cost  matrix,  or  when  we  are  referring  to  the 
original  cost  matrix  in  (l)  we  will  ignore  the  subscript  and  write  Z (a)  instead 
of  Zc(a)  . 

If  a  is  the  assignment  in  (2) 

rt  =  {(jj»  ')»  •••>  (jn,  n) } 

is  another  assignment  and  it  is  known  as  the  reflection  of  n  . 

An  assignment  b  is  called  a  symmetric  assignment  (SA  in  abbreviation)  if 
b  =  b,  i.e.,  whenever  it  contains  a  cell  (i,  j)  it  should  also  contain  the  cell 
( j ,  i)  .  Let  K  denote  the  set  of  all  SA.  The  symmetric  assignment  problem  is 
the  problem  of  finding  the  minimal  cost  (optimal)  SA. 

The  special  case  of  the  symmetric  assignment  problem  in  which  the  diagonal 
cells  bear  infinite  cost  is  equivalent  to  the  optimal  matching  problem  in  graph 
theory  T5].  The  famous  Chinese  postman  puzzle  can  be  formulated  as  a  matching 
problem  and  hence  as  a  special  case  of  the  symmetric  assignment  problem  (pointed 
out  by  Professor  D.  Gale). 


Terminology 


D iagona 1  Cel  1 :  Any  cell  along  the  principle  diagonal,  i.e.,  a  cell  of  the  form 
(i,  i )  is  cal  led  a  diagonal  cell. 

Node .  A  node  is  a  subset  of  K  of  the  following  form 

node  N  =  fall  SA  which  conta in  ( i ,  ,  j  , )  ( j ,  ,  i , )  ...  ( i  ,  j  )  ( i  ,  i  )  and  wh i ch 

1111  r  r  r  r 

do  not  contain  (m( ,  p 1 )  (p , .  m,  )  •••  (m^ ,  Ps)(p^,  m  )}  .  (3) 

Then  the  cells  ( i  j  •  j j )  ( j j  .  ij)  (*r>  jr)(j  >  'r)  are  said  to  be  the  cells 

specified  to  be  contained  in  the  node,  and  the  cells  (m^  ,  p ^ ) ( p ^ ,  m^ )  ...  (m^,  p  ) 

(pc ,  m^)  are  said  to  be  the  cells  specified  to  be  excluded  crom  the  node.  The 
letters  M,N  will  be  used  to  denote  the  nodes. 

For  simplicity  we  can  write  down  the  node  N  in  equation  (?)  as 

N  =■  f  ( '  1  .  jjMjj.  ',)  •••  ( '  r  -  JrHjr.  'r);(m,  .  P,)(P1  ~>Y)  • • • (m..  P?)  (p^TnT)} 

the  bar  above  a  cell  indicnt'ng  that  it  is  specified  to  be  excluded  from  the  node. 

n  admissible  unspecified  cell  at  node  N  is  any  cell  which  is  unspecified  at  node 
N  and  which  does  not  lie  in  a  row  or  column  of  a  cell  specified  to  be  in  that  node. 

Branching  from  node  N  with  the  admissible  unspecified  cell  (i,  j).  Suppose 
node  N  is  given  b'(  (3)  and  let  (i,  j)  be  an  admissible  unspecified  cell  at  N  . 
Then  ne i the  r  i  nor  j  equa I  any  of  ij,  ...,  i  ^ ,  j  ^ ,  ...,  j  .  Then  we  can 
part  1 1 ion  N  as 


-  '  -*»a 


k 


N  =  Nj  U  N^,  Np  disjoint,  where 

N,  =  j,)(j|.  '  J )  •••  ('r.  jr)(jr>  j)(j.  • )  ;  ("»|  .  Pj)(Pj,  m  )...(ms,  pg)  (p"  •  ms)) 

N2  ~  ( (  '  l  ’  j | ) ( i i  *  ' t )**"('  r  ’  j  ) ( j  r »  ir)(m]>  P| ) (p, .  PS)(P$»  nT)(T  j )  ( j  .  ')}• 

This  operation  of  partitioning  into  two  disjoint  subsets  is  known  as  branching  from 
node  N  with  the  admissible  unspecified  cell  (i,  j)  .  The  nodes  and  will 

be  called  the  branches  emanating  from  node  N  . 

Reducing  the  matrix  C  ,  method  I:  This  operation  proceeds  as  follows 

(1)  Subtract  from  each  element  the  minimal  element  in  its  row. 

(2)  In  the  resulting  matrix,  subtract  from  each  e’ement  the  minimal  element  in  its 

column.  Suppose  this  leads  to  a  matrix  Cj  ,  then  all  the  elements  in  Cj 

are  nonnegative  and  each  row  and  column  of  C  contains  at  least  one  zero  element. 

r  +  rT  1 

I  v  j  T  j  | 

(3)  If  Cj  is  not  symmetric,  let  Cj  =  - ^ -  .  Then  C(  is  a  symmetric  matrix 

with  all  nonnegative  elements. 

(L)  Suppose  contains  some  rows  with  no  zero  element  at  all.  Pick  one  of  them, 

say  row  i^  .  By  symmetry  column  i^  also  has  no  zero  element. 

(5)  Subtract  the  minimum  element  in  row  i^  from  each  element  in  that  row. 

(6)  In  the  resulting  matrix  subtract  the  minimal  element  in  column  i^  from  each 
element  in  that  column. 

Since  the  operation  in  (5)  alters  only  the  diagonal  element  ( i  ^ ,  i^)  in 
column  ij  ,  at  the  end  of  step  (6),  either  a  pair  of  symmetrically  placed  zero 
cells,  say  (Ij,  j)  and  ( j ,  i j )  are  created  in  row  ij  and  column  i^  respectively 
or  the  diagonal  cell  (ij,  ij)  becomes  a  zero  cell.  Suppose  the  result  is  Cj  . 
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If  is  not  symmetric  let  = 


c„  +  cl 

2  2 


By  symmetry  contains  all  the  zero 


cells  of  Cj  and  at  least  one  additional  zero  cell  in  each  of  row  i^  and  column 

I 

ij  .  If  C2  has  some  rows  and  columns  without  any  zero  entries  repeat  steps 

I  I 

(*0 .  (5) ,  (6)  with  C2  replacing  . 

Repeat  this  process  until  finally  a  symmetric  matrix  C  is  obtained  which 

K 

consists  of  all  nonnegative  elements  such  that  each  row  and  column  of  C  contains 

R 

at  least  one  zero  element. 

CR  is  known  as  the  reduced  matrix  obtained  from  C  .  The  sum  of  all  the 
numbers  subtracted  from  the  rows  and  columns  during  the  various  steps  is  known  as  the 
reduction  of  the  matrix  C  . 

Reducing  the  matrix  C  ,  method  2:  Using  this  method  the  operation  of  reducing 
C  proceeds  as  follows. 

(1)  Using  the  cost  matrix  C  find  the  optimal  assignment  by  the  Hungarian  method  [j  ]. 

The  Hungarian  method  transforms  C  by  adding  constants  to  its  rows  and  columns 
until  finally  a  matrix  Cj  is  obtained  which  consists  of  all  nonnegative  elements, 
with  at  least  one  zero  in  each  row  and  column.  The  optimal  assignment  is  contained 
among  the  zero  cells  of  . 

T 

•  C,  -C, 

(2)  If  Cj  is  not  symmetric,  then  repeat  step  (1)  with  ^ -  in  place  of  C  . 

Repeat  this  process  as  many  times  as  necessary  until  the  final  transformed 
matrix  with  nonegative  elements  and  containing  at  least  one  zero  in  each  row 


and  column,  is  symmetric. 
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The  sum  of  the  costs  of  all  the  optimal  assignments  as  they  are  obtained  in 
the  successive  steps,  in  known  as  the  reduction  of  the  matrix  C  .  The  final 
transformed  matrix  is  known  as  the  reduced  matrix  obtained  from  C  .  It 

can  be  seen  that  the  operation  of  reducing  C  by  method  2  might  involve  some  more 
work  than  that  by  method  1.  However,  the  reduction  of  C  obtained  by  method  2  is 
likely  to  be  larger  than  that  obtained  under  method  1  and  this  helps  in  improving 
the  efficiency  of  the  algorithm. 

It  has  been  proved  in  theorem  6  that  the  repetition  of  steps  (l)  and  (2) 
of  method  2  have  to  be  carried  out  only  a  small  number  of  times  before  obtaining 
.  Hence  is  obtained  from  C  in  a  small  number  of  steps. 

The  remaining  cost  matrix  at  node  N  :  Suppose  node  N  is  given  by  equation 
(3).  Then  the  matrix  obtained  by  striking  off  the  rows  and  columns  i^,  ....  i  j 

j  r  from  C  and  replacing  the  cost  elements  in  the  cells  (m^ ,  p^Hp^,  m^ )  .. 
(m  .  p$ )  ( Ps  ms)  by  infinity  (or  a  very  very  large  positive  number)  is  known  as  the 
remaining  cost  matrix  of  node  N  .  The  reduced  matrix  obtained  from  this  matrix  is 
known  as  the  reduced  remaining  cost  matrix  at  node  N  and  is  denoted  by  CN  R  . 

The  remaining  cost  matrix  of  K  is  C  itself,  and  the  reduced  remaining  cost 
matrix  of  K  is  . 

The  evaluation  of  an  admissible  unspecified  cell  at  a  node.  Let  ( i ,  j )  be 

an  admissible  unspecified  cell  at  a  node  N  .  Then  its  evaluation  at  node  N  is 

defined  to  be  9  (i,  j)  where 
N 

0N('.  ')  =  Sum  of  the  minimal  elements  in  row  i  and  column  i  of  R 
excluding  the  element  in  (i,  i) 


after 
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eN  (i.j) 


e  Sum  of  the  minima!  elements  in  rows  i  and  j  and  columns  i  and 
J  'n  C.  -  after  excluding  the  (i,  j)th  and  (j  ,  i)th  elements; 
if  these  minima  occur  at  distinct  places 
=  the  diagonal  element  in  (i,  i)  plus  the  sum  of  the  minimal  elements 
in  row  j  and  column  j  in  ^  after  excluding  the  (i,  j ) t h  and 
(j ,  i )  th  elements;  if  minimum  in  row  i  end  column  i  after  the 
exclusion  occurs  at  the  diagonal  cell  (i,  i)  and  the  diagonal 
element  at  (j  ,  j)  is  not  the  minimal  in  row  j  and  column  j  after 
the  exclusion 

=  the  sum  of  the  diagonal  elements  (i,  i)  and  ( j  ,  j )  in  ^  if 
those  elements  are  the  minimal  ones  in  rows  i,  j  and  columns  i,  j 
after  excluding  the  cells  (i,  j)  and  ( j ,  t )  . 


Since  the  reduced  matrix  contains  all  nonnegative  elements,  ©N('*  j )  >  0 
for  all  admissible  unspecified  cells.  Also  since  'he  reduced  matrix  has  at  least 
one  zero  element  in  each  row  and  column,  0^(i,  j)  =  0  unless  (i,  j)  is  a  zero 


cell  in  C 


N.R  ' 


The  mathematical  theory 

Lemma  1:  If  n  is  even,  the  total  number  of  SA  which  do  not  contain  any  diagonal 

ii  •  i  /  »n/2 
ce l 1 s  i s  n ! /  2 


Proof:  An  SA  which  does  not  contain  any  diagonal  cell  consists  of  pairs  of  cells 

like  (I,  j ) (J  ,  I),  1  j  .  When  n  is  even  the  total  number  of  ways  In  which  n 

n/2 

objects  can  be  paired  In  this  manner  is  n  1/  2  which  is  therefore  equal  to  the 
total  number  of  SA  without  any  diagonal  cells. 


( *  ■'M'n 


Theorem  1 :  The  total  number  of  possible  SA  is 


£  (n  -  r)  1/ ,  if  n  is  odd 


r=l 
r  odd 


i  M  (n 


r)  1/2 


,  /(,(n  -  r )/2  .  , 

,  i  f  n  is  even. 


r=0 
r  even 
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Proof:  This  follows  by  using  lemma  1  and  the  fact  that  if  n  is  odd,  any  SA 
must  contain  an  odd  number  of  diagonal  cells  and  if  n  is  even,  any  SA  contains 
an  even  number  of  diagonal  cells. 

Theorem  2:  Let  Cj  be  the  matrix  obtained  by  adding  a  constant  t  to  each  element 
is  a  row  (or  column)  of  C  .  If  b  is  an  optimal  SA  w.r.t.  cost  matrix  C  then 
it  is  also  an  optimal  SA  w.r.t.  cost  matrix  and  vice  versa. 

Proof:  This  is  actually  a  restatement  of  a  similar  theorem  for  the  general  assign¬ 
ment  prob’em  [1],  This  follows  easily  because 

Zc  (a)  *  Zj. ( cl)  +  l  ,  for  all  assignments  a  , 

and  t  is  a  constant. 

Theorem  3:  If  C  is  not  symmetric  the  optimal  SA  w.r.t.  cost  matrix  C  is  also 
optimal  w.r.t.  the  cost  matrix  C1  =  (C  +  C^)/2  . 

Proof:  This  follows  from  the  fact  that 

Z^(b)  =  Zj. ,  (b)  ,  for  all  beK 


V 
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Theorem  4:  The  reduction  of  C  is  a  lower  bound  on  the  cost  of  any  SA. 

Proof:  If  beK  then  by  theorems  2  and  3 

Zr(b)  =  reduction  of  C  +  Z  (b)  . 

L  CR 

But  since  each  element  in  CD  is  >0  we  have  Z  (b)  >  0  .  QED 

R  -  CR  — 

In  general,  corresponding  to  any  node  N  let 

LB(N)  =  original  cost  elements  of  all  the  cells  specified  to  be  contained  in  the 

node  +  reduction  of  the  remaining  cost  matrix  at  node  N  .  (4) 

Then  from  theorem  4  we  see  that  LB ( N)  is  a  lower  bound  on  the  cost  of  any  SA  in 
N  .  We  have  the  following  theorem  on  the  lower  bounds  of  the  branches  after 
branching  from  a  node. 

Theorem  5:  Let  (i,  j)  be  an  admissible  unspecified  cell  at  node  N  ,  which  is  a 

zero  cell  in  C..  D  .  Let  N.  ,  N.  be  the  branches  when  N  is  partitioned  w.r.t 
N ,  K  I  Z 

(  i  .  j  )  • 

=  {beN  :  (i ,  j)eb,  ( j ,  i ) eb } 

N2=  {beN  :  (i ,  j)^b,  (j ,  i)^b}  . 

Then  (i)  LB(N2)  >  LB(N)  +  0N{i,  j) 

(ii)  Let  C  be  the  matrix  obtained  by  striking  off  the  rows  i,  j  and  columns 
N1 

j,  i  from  CN  R  and  let  be  the  matrix  obtained  by  replacing  the  entries  in 

the  cells  ( i  ,  j)  and  (j ,  i)  of  CN  R 
number).  Then 


by  infinity  (or  a  very  large  positive 
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LB(N.)  =  LB(N)  +  reduction  of  Ck, 

1  N, 

LB(N-)  =  LB(N)  +  reduction  of  CN 

2 

P roo f :  (i)  This  follows  from  the  fact  that  by  definition  beN^  ^  beN  and 

( i  .  j )  4b  ,  (j  ,  i)  ^b  . 

(ii)  These  follow  the  fact  that  (i,  j)  which  is  admissible  unspecified  at 

N  i s  a  zero  cel  I  in  C..  . 

N  ,  K 

Theorem  6:  In  reducing  any  matrix  C  by  method  2,  step  2  need  not  be  applied  more 
than  twice. 


Proof:  Suppose  C.  is  the  transformed  matrix  obtained  when  step  1 

•  C1 

C  .  If  Cj  is  symmetric  we  are  done.  Otherwise  we  find  C(  =  - 

a  symmetric  matrix  wirh  nonnegative  elements. 


i s  app lied  to 

+  cT 


1 


which  is 


Let  a  be  an  optimal  assignment  w.r.t.  cost  matrix  C(  .  Let  be  the 

I 

transformed  matrix  obtained  when  step  (l)  of  method  2  is  applied  to  Cf  .  Then 
for  any  assignment  b 


But 


Hence 


Z  ,  (b)  =  Z  ,  id)  +  Z-  (b)  (by  the  Hungarian  method) 


Z  ,  (a)  =  Z  ,  ia)  because  C  is  symmetric  . 


Zr  (a)  =  Zr  (a)  =  0  . 
2  L2 


Since  has  all  nonnegative  elements  this  implies  that  the  entries 


in  the  matrix  in  the  cells  of  a  and  a  are  all  zero.  This 


1 1 


c  + 

I  ^2  C2 

implies  that  in  =  - - -  there  exists  at  least  one  zero  cell  in  each  row  and 

I 

column  (corresponding  to  the  cells  of  a  and  a)  and  hence  0^=0^  . 

Theorem  7:  Let  N  be  any  node.  Let  a  be  the  minimal  cost  assignment  in  N  . 

(i)  If  N  contains  only  a  single  assignment,  then  it  must  be  a  and 
Z (a)  =  LB(N)  . 

(ii)  If  method  2  is  employed  for  reducing  matrices  then  Z(a)  =  LB(N)  . 

Proof:  (i)  This  follows  from  the  way  LB(N)  is  defined  in  equation  (4). 

(ii)  Z(b)  =  LB (N)  +  Zr  (b) 

LN,R 

for  all  assignments  b  in  N  ,  by  theorem  5  and 

Z _  ( a)  -  0  (by  Hungarian  method) 

LN,R 

Therefore  Z(a)  =  LB(N)  . 

Algorithm: 

Stage  1  :  Find  LB(K)  .  If  this  is  found  by  using  method  2  for  reducing  C  ,  the 
algorithm  terminates  if  there  exists  an  optimal  assignment  which  is  symmetric. 
Otherwise  branch  from  K  . 

General  stage  m:  At  this  stage  K  has  been  partitioned  into  several  nodes.  Any 
node  which  has  not  been  branched  is  called  a  terminal  node  at  this  stage.  In  the 
course  of  the  algorithm  a  lower  bound  on  the  cost  of  the  minimal  SA  in  each  node 
(corresponding  to  equation  (4))  has  been  obtained.  A  terminal  node  which  has 


/ 
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the  least  value  of  the  lower  bound  among  all  the  terminal  nodes  is  known  as  a 
minimal  terminal  node  at  this  stage. 

Optimality  criterion:  The  algorithm  is  terminated  whenever  there  exists  a  minimal 
terminal  node  which  contains  only  one  SA  or  when  there  exists  a  minimal  terminal 
node  in  which  the  optimal  assignment  is  an  SA  (if  method  2  is  used  for  reducing 
mat  r i ces) . 

If  the  optimality  criterion  is  not  satisfied  at  this  stage,  then 

(i)  find  out  the  minimal  terminal  node  with  the  least  cardinality  (i.e.,  which 
has  the  maximum  number  of  cells  specified  to  be  contained  in  it,  among  all  the 
minimal  terminal  nodes  at  this  stage).  Suppose  it  is  node  N 

(ii)  find  (i,  j)  ,  an  admissible  unspecified  cell  at  node  N  such  that 

(a)  (i,  j)  is  a  zero  cell  in  CN  R  and 

(b)  0  (i,  j)  is  maximum  among  all  the  cells 
satisfying  (a)  . 

(iii)  branch  from  N  w.r.t.  (i,  j)  and  find  the  lower  bounds  of  both  the 
branches  using  theorem  5> 

( iv)  go  to  stage  (m  +  1 )  . 

These  computations  are  repeated  until  the  optimality  criterion  is  satisfied. 

The  optimal  assignment  (or  the  only  assignment)  in  the  minimal  terminal  node  at 

the  final  stage,  which  is  an  SA,  is  an  optimal  SA. 
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Justification  for  the  algorithm:  At  each  stage  the  terminal  nodes  are  mutually 
disjoint  and  their  union  is  K  .  So  the  optimal  SA  has  to  lie  in  one  of  the  terminal 
nodes  and  it  is  more  likely  to  lie  in  a  minimal  terminal  node  than  in  any  other. 

This  is  the  reason  for  branching  from  a  minimal  terminal  node  at  each  stage. 

Also  if  the  optimal  assignment  (or  the  only  assignment)  in  a  minimal  terminal 
node  is  an  SA  then  that  SA  must  be  an  optimal  SA  by  theorem  7-  This  proves  the 
validity  of  the  optimality  criterion. 

It  can  be  seen  that  the  efficiency  of  the  algorithm  improves  if  we  could 
guarantee  that  we  reach  rapidly  a  minimal  terminal  node  which  contains  a  single  SA. 
When  a  node  N  is  branched  as  in  theorem  5  it  can  be  seen  that  the  cardinality  of 

is  far  greater  than  that  of  .  Hence  we  choose  the  cell  (i,  j)  for  branching 

N  ,  in  such  a  way  that  the  lower  bound  of  is  made  as  large  as  possible.  By 

( i )  of  theorem  5,  the  9^(i,  j)  can  be  used  to  achieve  this  propose  and  this  is 
what  we  used  in  step  (ii)  of  the  general  stage  of  the  algorithm.  The  reason  for 
branching  from  the  minimal  terminal  node  with  least  cardinality  is  also  to  help 
in  reaching  a  minimal  terminal  node  with  a  single  SA  rapidly. 

Finally  method  2  for  reducing  matrices  helps  in  obtaining  a  more  precise  lower 
bound  on  the  cost  of  the  minimal  SA  in  any  node.  And  it  helps  in  checking  whether 
the  optimality  criterion  is  satisfied  at  a  much  earlier  stage.  Thus  eventhough 
method  2  might  involve  more  computatinal  work  at  each  node,  it  may  require  much 
less  branching  before  the  optimal  SA  is  found.  Thus  the  use  of  method  2  for  reducing 
matrices  might  lead  to  a  more  efficient  algorithm,  especially  when  n  is  large. 
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A  numerical  example: 
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Method  1  has  been  used 
to  reduce  matrices. 

Lower  bound  corresponding 
to  each  node  is  indicated 
by  its  side. 
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The  Chinese  postman  problem: 

This  is  the  problem  of  finding  the  minimal  distance  route  passing  through  each 
edge  of  a  connected  undirected  graph  and  then  returning  to  the  origin.  Hence  such 
a  route  is  also  called  an  Edge  Covering  Tour  (E.C.T.  in  abbreviation).  It  is 

assumed  that  the  distance  associated  with  each  edge  of  the  graph  is  nonnegative. 

We  denote  the  graph  by  G  .  We  shall  call  the  nodes  of  the  graph  as  vertices  to 
avoid  confusion  with  the  nodes  defined  in  the  branch  and  bound  algorithm.  A  vertex 
of  the  graph  at  which  an  odd  number  of  edges  are  incident  is  known  as  an  odd  vertex. 
If  there  are  no  odd  vertices  then  there  exists  an  E.C.T.  in  which  each  edge  of  the 
graph  is  travelled  exactly  once.  Such  an  E.C.T.  is  known  as  an  Euler  cycle. 
Obviously  if  an  Euler  cycle  exists  it  gives  the  minimal  distance  E.C.T. 

If  the  graph  contains  some  odd  vertices  then  their  total  number  is  even.  Let 
them  be  numbered  1,  2,  ...,  2n  . 


Lemma  2 :  There  exists  a  minimal  E.C.T.  in  which  no  edge  of  the  graph  G  will 
be  travelled  more  than  twice. 


Proof:  Let  t  denote  the  minimal  E.C.T.  Obtain  a  new  graph  G  by  drawing 
each  edge  of  G  as  many  times  as  it  is  travelled  in  t  .  Then  t  is  an 

I 

Euler  cycle  of  G  . 


Hence  each  odd  node  of 

number  of  repeated  edges  In 
number  of  repeated  edges  In 

Now  obtain  a  new  graph 

I 

odd  number  of  times  in  G 


G  has  become  an  even  node  in  G  .  Therefore,  the 

I 

G  Incident  at  each  odd  node  in  G  is  odd,  and  the 

I 

G  incident  at  each  even  node  In  G  is  even. 


1 1 

G  by  drawing  each  edge  of  G  once  if  it  occurs  an 

I 

and  twice  if  it  occurs  an  even  number  of  time  in  G 
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I  1 

By  the  above  property  it  is  clear  that  all  nodes  are  even  nodes  in  G 

I  I 

So  '  an  Euler  cycle  t  which  is  an  E.C.T.  of  G  ,  and  the  total  distance 

ii  ii 

travelled  in  t  is  less  than  or  equal  to  that  in  t  .  Also  in  t  each  edge  of 
G  is  travelled  either  once  or  twice. 

Lemma  3:  There  exists  an  optimal  E.C.T.  ,  tg  ,  on  the  graph  G  with  the  following 

property.  No  edge  of  G  is  travelled  more  than  twice.  Let  T\  be  the  set  of 

0 

edges  of  G  which  are  travelled  twice  in  tg  .  Then  the  set  of  odd  vertices  of 

G  can  be  partitioned  into  pairs  (ijp  •••  ('n|>  'n2^  suc^  that  the  edges 

in  f  form  n  different  edge  disjoint  paths  connecting  each  pair  of  odd 
t0 

vertices  above. 


Proof:  Let  G.  be  the  graph  obtained  by  adding  the  edges  of  to  G  .  Then 

t0  ro 

tn  is  an  Euler  cycle  in  G  .  Hence  if  i  is  an  odd  vertex  in  G  ,  then  V 
U  *0  C0 
must  contain  a  path  from  i  to  some  other  odd  node  j  .  Suppose  along  this  path 

in  r  there  lie  two  other  odd  vertices  r,  s  .  Adding  all  the  edges  of  this 
t0 

path  to  G  leaves  r,  s  odd  still.  Hence  by  the  assumptions  in  the  hypothesis 

T  must  consist  of  another  path  between  r  and  s  .  But  then  the  E.C.T., 
t0 


tj  ,  obtained  by  deleting  all  the  edges  along  both  these  paths 
r  and  s  is  likely  to  be  better  than  tg  . 


i  n 


between 


Hence  we  can  assume  that  there  exists  an  optimal  E.C.T.,  tg  ,  with  the 

property  mentioned  in  the  lemma.  Also  in  any  of  the  paths  in  T  between  a  pair 

*0 

of  odd  nodes,  not  more  than  one  other  odd  node  can  lie. 


Lemma  4:  Let  tg  he  an  optimal  E.C.T.  with  the  property  mentioned  in  lemma  3* 


let  i,2) 


(inr  'n2 


)  be  a  corresponding  pairing  of  odd  vertices.  Then 


the  paths  in  1^  between  any  pair  ('r|«  i^)  must  be  a  shortest  route  between 
that  pair  (i  ,  i  )  on  G  . 
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P roo f ;  The  result  follows  obviously  because  if  the  path  in  T  between  the  pair 

l0 

of  odd  vertices  were  not  a  shortest  route  in  G  between  that  pair,  then  a  better 

E.C.T.  can  be  obtained  by  deleting  the  existing  path  in  T  and  replacing  it 

t0 

by  the  shortest  route. 

Lemma  2,  k  together  imply  that  the  optimal  E.C.T.  corresponds  to  a  minimal 
cost  SA  with  the  matrix  of  shortest  route  distances  between  odd  vertices  as  the 
cost  matrix.  [Pointed  out  by  Professor  0.  Gale.] 

Algorithm  for  the  Chinese  postman  puzzle: 


Let  c..  =  o'  ,  a  very  large  positive  number 


i  i 


c.  .  = 
•J 

i*j 


total  distance  of  the  shortest  route  on 

graph  G  from  odd  vertex  i  to  odd  vertex  j  , 

i ,  j  “  1 ,  2 ,  ....  2n  . 


C  =  ( c ..)  . 

Find  the  minimal  SA  corresponding  to  the  cost  matrix  C  .  If  the  minimal  SA  is 


(i,.  J  i )  (j,’  V  V  (jn’  'n) 


then  obtain  a  new  graph  by  duplicating  all  the  edges  along  the  shortest  routes  of 
( i  1  ,  j^)  ( i  ^  ,  j‘2)  ...  ( i  ^  ,  j  )  respectively.  The  new  graph  obtained  has  an  Euler 
cycle,  which  gives  the  optimal  E.C.T.  on  G  . 


The  proof  of  the  algorithm  follows  from  lemmas  2,  h. 
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